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Abstract 



A family of exact relativistic stellar models is described. The family generalizes Buchdahl's n = 1 
polytropic solution. The matter content is a perfect fluid and, excluding Buchdahl's original model, it 
behaves as a liquid at low pressures in the sense that the energy density is non-zero in the zero pressure 
limit. The equation of state has two free parameters, a scaling and a stiffness parameter. Depending on 
the value of the stiffness parameter the fluid behaviour can be divided in four different types. Physical 
quantities such as masses, radii and surface redshifts as well as density and pressure profiles are calculated 
and displayed graphically. Leaving the details to a later publication, it is noted that one of the equation 
0^ . of state types can quite accurately approximate the equation of state of real cold matter in the outer 

regions of neutron stars. Finally, it is observed that the given equation of state does not admit models 
^-l with a conical singularity at the center. 

U 

1 Introduction 



Among the new exact relativistic stellar models recently found in [jlj, 0] , two families are of particular physical 
interest. They both generalize gaseous models found by Buchdahl |gf Q. His models have equations of state 
which reduce to a Newtonian polytropic gas (of index n = 1 and n — 5 respectively) at low pressures. We will 
refer to the generalized solutions as the generalized Buchdahl n = 1 polytrope (GB1) and generalized Buchdahl 
n = 5 polytrope (GB5) models respectively. By contrast to Buchdal's original solutions the generalized models 
have equations of state which are liquid at low pressures in the sense that the density is non-zero at zero 
pressure, p s > Oj^] Their equations of state each have two free parameters. By properly adjusting those 
parameters, the models can describe the interior of neutron stars with surprising accuracy In this paper 
we describe the GB1 family 0. The matter in the model is a perfect fluid obeying a certain barotropic 
equation of state p = p(p). The history of static spherically symmetric (SSS) models with matter goes back 
to Schwarzschild's discovery already in 1916 of his interior solution with incompressible matter |5fl. Although 
a number of exact SSS models have been found since then (see e.g. Tolman's classic paper ||) there is an 
important distinction which should be made between solutions such as Schwarzschild's interior solution and 
solutions of Tolman's type. Whereas the former kind is a family of solutions describing a stellar sequence 
with varying masses for the given equation of state (p = constant in that case) the latter kind gives only 
a single solution with a fixed mass once the equation of state is specified. The underlying reason is that 
Schwarzschild's interior model is the general solution of the Einstein field equations for the given equation 

1- We use a subindex "s" throughout the paper to denote the surface values (where p = 0) of various physical quantities. 
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of state. By general solution we mean a solution with the full number of integration constants. For any 
given equation of state there is a family of SSS models parametrized by precisely two physically relevant 
integration constants. Usually one of these constants is determined by requiring that the stellar model is 
regular at the center. The single remaining nontrivial integration constant then determines the total mass. 

The general solutions correspond to nontrivial symmetries of the field equations. Such symmetries only 
exist for certain special equations of state [jj]. Tolman type solutions on the other hand are particular 
solutions which do not contain the full number of integration constants. Such special solutions have been 
referred to as submanifold solutions since they correspond to submanifolds in the phase space which are left 
invariant by the field equations |Q. Obviously the general solutions are more useful than the submanifold 
solutions. However, until recently Schwarzschild's interior solution was the only general SSS solution which 
was considered to be of physical interest. 

At this point we should make clear what we mean by a physically interesting solution. For the purposes 
of this paper we consider a solution to be physical if the energy density is non-negative, p > (gravity is 
attractive), and that the matter is locally mechanically stable so that p > and dp/ dp > 0. Incompressible 
matter is considered to be included in this category. Schwarzschild's interior solution was the first matter 
filled general physical SSS solution of the Einstein equations. It took almost half a century before the 
next general physical SSS model, the n — 5 polytrope solution, was discovered by Buchdahl in 1964 0. 
As in the case of a Newtonian n — 5 polytrope this solution has an infinite radius although the mass is 
finite. This model was not considered to be of any significant physical interest. A few years later Buchdahl 
found the n = 1 polytrope solution ||. Although this model has a finite radius it has not been widely 
discussed in the literature. The next development came recently when the generalizations of Buchdahl's 
solutions were discovered by Simon § and Rosquist @. The GB1 and GB5 solutions are in fact families of 
solutions with equations of state depending on two continuous parameters. In the present paper we discuss 
the GB1 family in some detail. Because of the two parameters the structure of this family is quite rich with 
qualitative behaviour varying in different parts of parameter space. In particular, by an appropriate choice 
of the equation of state parameters, the GB1 family can be fitted quite accurately to the Harrison- Wheeler 
(HW) equation of state [|| for cold matter expected to be valid in the outer parts of a neutron star. Later 
refinements of the HW equation of state are essentially corrections in the regime above nuclear densities. 
The fitting will be performed in a separate paper ||. It turns out that the best fit GB1 model is very close 
to the original Buchdahl solution rendering that model more interesting than previously thought. As for the 
GB5 family, its basic properties are planned to be discussed in a forthcoming paper [jl0| (see also pi) ) . 

One reason that the new static models are interesting is that they can serve as starting points to find an 
interior solution for the Kerr metric. Assuming that such an interior general physical solution exists, and 
that it is valid for a slowly rotating Kerr exterior, there must be an SSS interior solution in the limit of zero 
rotation. Furthermore, this static limit will also be a general physical solution because of the arbitrary mass 
parameter which is present in the Kerr metric. Of course this heuristic argument excludes interiors with 
singular matter distributions such as the Negebauer-Meinel disk Jl2| . Referring to the extensive classification 
of general SSS models in Q, the GB1 and GB5 solutions are in fact the only good candidates for the static 
limit of a possible interior Kerr solution. On the other hand it is quite conceivable that there are rotating 
generalizations of the GB1 and GB5 models (as of any fluid SSS model) which have an exterior vacuum 
field which is different from Kerr. One might try to spin up the GB1 model to yield a rotating perfect fluid 
solution which may or may not have an exterior gravitational field which coincides with the Kerr solution. 
Attempts in this direction have been made by trying to extend the use of the Newman- Janis trick |ll| to 
impart rotation to fluid SSS models |lj, [ll| but without success so far. However, the GB1 or GB5 models 
so far have not been considered in this context. 

Exact SSS solutions may also serve as benchmarks to test the numerical codes for integrating the equations 
of stellar structure. Such testing was performed for example in |l6| using the exact solution for the n = 1 
Newtonian polytrope. The existence of exact relativistic solutions with physically reasonable equations of 
state make it possible to test such codes also in the relativistic domain. 

Another important use for exact solutions is to investigate in more detail how the maximum mass depends 
on properties of the equation of state. Only general physical solutions qualify as objects to study in this 
context since it is necessary to vary the mass while keeping the equation of state fixed. However, in the past 
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there was no exact model available which was capable of exhibiting a mass maximum. The reason is that the 
mass maximum requires a rather soft equation of state by contrast to the hardest conceivable matter, the 
incompressible perfect fluid, which makes up the matter source for the interior Schwarzschild solution. The 
GB1 models which are the subject of this paper also contain matter which is too hard for a mass maximum 
to appear. However, as will be shown in |10[ ], the GB5 models have a softer equation of state and do exhibit 
a mass maximum. 

In addition to the topics discussed above the GB1 models have also turned out to be interesting from the 
point of view of understanding the space of solutions of the SSS Einstein equations. Taking the case of an 
incompressible fluid as an example, one of the two physically relevant integration constants is determined 
by the requirement that the metric should be regular at the center. This is sometimes referred to as the 
condition of elementary flatness. The center is defined as that value of the radial parameter for which the 
area of the 2-sphere of constant radius tends to zero. However, as described in section || the GB1 models 
have a different behaviour in this respect. It turns out that a GB1 model can only be continued to the 
center for one specific value of one of the integration constants. The condition of elementary flatness is then 
automatically satisfied. The global behaviour of singular GB1 models therefore seem to be different from 
that of singular incompressible models. A discussion of possible singularity types for SSS models was given 
in |l7j but we have not made an attempt to see how the GB1 models would fit in that scheme. 

2 The model 

The usual starting point when studying SSS models is to write down the metric in Schwarzschild coordinates 

ds 2 = _ e 2u dt 2 + e 2X dj ,2 + r 2 dn 2 ^ (1) 

where 

dfl 2 = d8 2 +sin 2 8d<j) 2 , (2) 

is the metric of the 2-sphere and v and A are functions of Schwarzschild's radial variable, r. However, this 
form of the metric is not general enough in the sense that it is not possible to express all exact solutions 
in closed form using Schwarzschild coordinates. It is therefore useful to write the SSS metric in terms of a 
general radial variable, R, 

ds 2 = -c 2l/ dt 2 + N 2 dR 2 + W 2 dn 2 , (3) 

where N and W are both functions of R. Schwarzschild's radial variable then becomes a derived quantity 
given by r — W(R). The form of the metric (|J) reflects the free gauge choice of the radial variable expressed 
by the freedom to choose the radial gauge function N. To write down the GB1 model we use the gauge 
N = e~ v . The SSS metric then becomes 

ds 2 = -Zdt 2 + Z^dR 2 + W 2 d(l 2 , (4) 

where W and Z := c 2v arc functions of the radial parameter R. Since N > 0, the radial variable R is 
increasing outwards from the center of the star. The Schwarzschild radial variable is given by r — W(R). 
The center of the staij^ is therefore located at R = R c defined by the condition W(R C ) = 0. The GB1 model 
is defined by § | 

Z = V/W , W = T + X , V = T - X , 

T= f(^-)- 1 (coshC)sinK(i?- J R_)] , for5>l , 

[(coshC)OR--R-) for 5 = 1 , U 

X = (cj + )- 1 (sinhC)sin[w + (i?-i? + )] , 



2 We use the subindex "c" to denote the center of the star. 

3 For convenience we have performed the rescalings t — > 4t, ds 2 — > 16ds 2 , T — * 4T and X — * 4X relative to the conventions 
of reference El . 
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where u>± = yj2na{8 ± 1) and C, > 0. The parameters £ and R± are integration constants while the 
parameters a and 8 characterize the equation of state. The latter can be parametrized as 

p = a{Z 2 - 26Z + 1) , p = a(-5Z 2 + 66Z - 1) , (6) 

where a > and 8 > 1. As discussed in models with (5 < 1 also exist but can only be used in the pressure 
region p > a(l — S 2 ). The range of <5 will be further restricted below. From the physical point of view, a is a 
scaling parameter. Models with different values of a can be transformed into each other by simply scaling the 
metric. The parameter a represents an overall scaling of length scales in the model. A model with given a 
has a natural length scale given by I — {na)^ 1 / 2 . However, scale invariant physical quantities such as surface 
redshift are unaffected by changes in a. As will be explained below the parameter 8 can be interpreted as 
a measure of the stiffness of the matter. We will refer to it as the stiffness parameter. The integration 
constants are not all physical. Only £ and the difference D := R + — i?_ carry physical information. This 
is because D is invariant under translations of R. Because of its nature as a metric coefficient Z must be 
strictly positive. Physical requirements on the equation of state further restrict the range of Z as discussed 
in the next section. 

Expressing the functions T and X in terms of r and Z we have 

T = ±r(l +Z), X = ±r(l - Z) . (7) 

Since the conditions r > and Z > must hold in the stellar interior it follows that T > 0. Also, as shown 
below in ([l2]), we have that Z < 1 which shows that X > 0. 



2.1 The equation of state 

In this paper we assume that the energy density is non- negative, p > 0, and that the matter is mechanically 
stable so that p > and dp /dp > 0. It follows immediately from (Js[) that the condition p > is satisfied 
precisely if 



Z < Z s = 6 - yjK 2 - 1 , (8) 
where Z s denotes the value of Z at the star surface defined by the condition p = 0. Also, the relation 

dp 5 — Z 



dp 5Z — 3(5 

together with Z > implies that the condition dp /dp > is equivalent to the relation 

36 



(9) 



„ < z < S . (10) 
5 



The speed of sound is given by u := y/dp/dp. It reaches the speed of light, u = c = 1, at Z* = 26/3. 
It is sometimes of interest to allow superluminal sound speeds for example to obtain a simple model such 
as Schwarzschild's interior solution. The limit in which the speed of sound becomes superluminal will be 
referred to as the causal limit and denoted by an asterisk subscript. The pressure at the causal limit is given 
by p* = a(l — 86 2 /9) and the corresponding energy density is p* = a(166 2 /9 — 1). At zero pressure the 
energy density is p s — 4a(l — 5 2 + 5yj6 2 — 1). The speed of sound at the surface is given by 



Ut = ; . (11) 

It follows from this expression that the zero pressure sound speed reaches the speed of light at 6 = 8* := 
3\/2/4 w 1.061 and becomes infinite in the limitQ 6 — ► <5oo := 5/V21 ~ 1.091. Thus models with 6 > 8* have 

4 We denote the limit of infinite sound speed by a subscripted infinity symbol. 
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a supcrluminal speed of sound for all physcially allowed values of the pressure while models with 8 < 8* are 
causal for p < p* and acausal or p > p* . 



The upper limit in must be strictly larger than the lower limit in (10). This gives the additional 
restriction 8 < 8oo and then 

3<5 



, < Z < 8 - y/8 2 - 1 < 1 . (12) 
o 

It can be checked that the condition p > is always satisfied if 8 > 1 and Z is in the range (|l2|). 

An important physical characteristic of the equation of state is the relativistic adiabatic index defined by 

It can be interpreted physically as a compressibility index. Lower 7-values correspond to softer matter. For 
the GB1 model the adiabatic index is given by 

7 = 4 ^~^ 2 . (14) 

1 (5Z -3S)(Z 2 -2SZ + 1) 1 ' 

For models with 8 > 1 the adiabatic index tends to infinity at the surface while 7 S = 2 for Buchdahl's original 
model (S = 1). The adiabatic index also tends to infinity in the limit Z — > := 3<5 / 5 corresponding to 
Vsound — > 00. At the causal limit the adiabatic index is given by 7* = 8<5 2 /(9 — 88 2 ). 

In practical calculations it is convenient to use £ :— ui-/oj + = y/(8 — l)/(8 + 1) as a stiffness parameter 
in place of 8. The value £ = (Buchdahl's original model) corresponds to the softest equation of state within 
the GB1 family. The range of £ is 

< £ < £00 := ^(5- \/21) w 0.2087 . (15) 

The causal limit is at £ = 3 — 2y2 ~ 0.1716. The parameter £ can be interpreted physically in terms 

of the speed of sound. Using (^) and (0) and defining fe : — p/p we find that 

{ =V3TTv5' (16a) 



where 



d:=d 1 -d 2 = 4u 2 (u 2 -2k- 5u 2 k) , (16b) 
dx := (fc + 1)(1 + 5m 2 ) 2 , (16c) 
d 2 := k + 1 + 2(9fc + b)u 2 + 3(15fc + 7)u 4 . (16d) 



In the limit fc — > this reduces to 

£= ; 2m ' = = ^ s 2 + 0«). (17) 

1 + 5u 2 + ^1 + 10u 2 + 21uj 

The inverse of this relation takes the simple form u 2 = where we have defined the parameter /i := 
1 — 5£ + £ 2 which takes values in the interval < p < 1. Thus for u s <C 1 we see that £ is essentially the 
sound speed squared. Assuming that the material in a cold star at low pressure consists of pure iron we can 
get a rough estimate of the value of £ for a realistic stellar model. The speed of sound in iron is 6 • 10 5 cm/s 
Jl8[ corresponding to u = 2 • 10~ 5 in geometrized units. Inserting this value in ( fl6| ) gives £ ~ 4 • 10~ 10 . This 
value is not sensitive to the value of p for moderate pressures up to and beyond p — 1 atm. In [|| we will 
determine £ more precisely by fitting the GB1 equation of state to the HW equation of state. This will also 
involve a determination of the scaling parameter a. 
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Figure 1: The two equation of state parameters xo and x\ are plotted against £ using logarithmic scales. Both xo 
and x\ tend to infinity in the limit £ — > ^ w 0.2087. They are both approximately unity at £ = 0.06. 



It is also convenient to replace Z by a new normalized parameter, x, which takes values in the interval 
[0,1] with x s = and Xoo = 1. This is accomplished by the transformation Z — Kx + Z s where K := 
Zoo - Z s = -(2/5)/V(l - £ 2 ) and Z s = (1 - 0/(1 + which implies 

p = &x(x + xo) , (18a) 

p = bb{-x 2 + 2x + x 1 ) , (18b) 

where 6 = aif 2 and 

x = mti- 1 , X! = 10£(1 - £) V 2 . (19) 

These constants satisfy the inequality < xq < x\ . They are plotted as functions of £ in figure [|. We see 
that xq and x\ are approximately equal up to £ ~ 0.06 corresponding to xo,xi < 1. In the expression for 



t he p ressure (18a) the term xq comes into play if xq > x. Similarly in the expression for the energy density 
(18b) the term x\ is important if x\ > x(2 — x). Since the functions x and x(2 — x) both take values in the 
interval [0, 1], this leads to a rough subdivision of the GB1 equations of state into the four types 

(I) x = x x = (£ = 0) , 

(II) 0< Vl «l (0<£< 0.001), 

(III) xo, xi - 1 (0.001 < ^ < 0.2) , 

(IV) x , xi > 1 (0.2 < ^ < ^oo « 0.2087). 

As discussed above, comparison with the actual equation of state of cold matter at low pressure demands 
£ ~ 4 • 10~ 10 corresponding to a type II equation of state. The four different types of equation of state are 
plotted in figures || and ||. The type I and II equations of state are practically indistinguishable in figure || 
but differ drastically in the behaviour of 7 in the zero pressure limit as illustrated in figure [5| 

In figure I the pressure and energy density are plotted as functions of x and normalized by the energy 
density at the causal limit which occurs at 

x-x - 5(1-6* + ^) (21) 
* 6(1-5^ + ^) ■ (21) 

The energy density and pressure in that limit are given by 

= §1^-2(7 + e 2 )(l + 7e 2 ) , p. = §^- 2 (l-34£ 2 +£ 4 ) . (22) 
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Figure 2: The pressure and the energy density are plotted for four values of £ representing the different types of 
equation of state as discussed in the text. For each of these values the upper curve corresponds to the normalized 
energy density while the lower curve gives the normalized pressure. The causal limit p — p, in this diagram is the 
line p=l- 



Since any given value of p corresponds dimensionally to a pressure pc 2 it is reasonable to normalize both 
the energy density and the pressure by . The dimensionless normalized pressure and energy density then 
become 

p :=vlP* = bx{x + x ) , ^ 
p := p/p* = 5b(—x 2 + 2x + xi) , 



where 

36(1 -5£ + <^ 2 
25(7 + £ 2 )(l + 7£ 2 

To see the full picture when plotting 7 we use the pressure p^ at which 7 tends to infinity as a reference 
point. Setting x = 1 in (18a) gives 

1 + 5£ + e 
1 - 5$ + e 



b — y ^il^ (24) 



Poo=b- - ^ ; ^ . (25) 



Then we have 



where 



p/Poc = bx(x + x ) , (26) 
b + e (27) 
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Figure 3: The adiabatic index for the n = 1 Buchdahl equation of state and some of its generalizations. The 
pressure is scaled by the pressure at which the adiabatic index becomes infinite. For £ > the adiabatic index 
diverges in the zero pressure limit. The indicated causal limit is the boundary between the regions with subluminal 
and superluminal sound velocity. It is seen that the causal region defined by subluminal sound speed is located in the 
lower left hand corner of the diagram with a narrow section stretching towards infinity along the 7-axis. Models with 
£ < = 3 — 2y/2 « 0.1716 start in the causal region at low pressures and enter the acausal region at the pressure p* 
while the models with £ > are acausal for all pressures. 



Writing the adiabatic index in terms of the normalized parameter gives 



2(5£ + Mz) 2 [M5-2:r) + 15fl 



For Buchdahl's polytrope this reduces to 



■<= 2 -m- 



At zero pressure we have 7 = 2 for this model. For the models with £ > the adiabatic index instead tends 
to infinity in that limit with the asymptotic form 7 ps /? s itgP -1 . 
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3 Physical characteristics of the stellar model 



In this section we compute some of the general physical characteristics of the GB1 family. As already stated 
the actual fitting to the HW equation of state will be done in the sequel paper Q . First we briefly discuss 
criteria which determine the degree to which the star may considered as relativistic. One such criterion is 
the surface redshift z s which in principle can be measured directly. In practice, however, this seems to be 
unrealistic for neutron stars. However, redshifts can be determined indirectly if the mass and radius of the 
star is known. The surface redshift is given by the formula z s = (1 — 2/3) -1 / 2 — 1 where := M/r s (see 
e.g. |^9|). A given value of the redshift corresponds to an equivalent velocity through the Doppler relation 
1 + z = v(l + v/c)/(l — v/c). For a Newtonian star, z s « 1. Defining a star to be relativistic if the surface 
redshift corresponds to a velocity greater than one percent of the speed of light this would then be the case 
if z > 0.01. In principle the redshift can attain values from the Newtonian z s = up to Buchdahl's limit 
[pp[ z s = 2 (corresponding to r s > jM) which is the maximum redshift if the density is assumed to be 
non-increasing in the outward direction, dp/ dr < 0. In practice, however, the surface redshift of a neutron 
star is expected to be in the range 0.1 < z s < 0.9 pi]] . Another criterion can be formulated by considering 
the deviation from Euclidean geometry. This deviation can be defined as the quantity a := 1 — r/ R pmp 
where i? pr op is the proper radial distance defined for a general SSS metric by i? pr0 p = J R NdR. The star 
can then be said to be relativistic if, say, a s > 0.01. 

Regular stellar models must satisfy the condition for elementary flatness at the center. This can be 
conveniently expressed in terms of the Schwarzschild radial gauge function which for the GB1 models is 
given by @] 

^Schw = Z-V^Wm)- 1 . (30) 
The elementary flatness condition then takes the form 

lim N Schw = 1 . (31) 

Using ( |30| ) it can be written as Z^ 2 (W,r) c = 1. From (||) we find (X tR ) c = cosh£, (Jjj) c = cosh£. This 
implies (W, H ) C = e<, {V, R ) C = e"< and Z c = (V/W) c = (V R ) C /(W M ) C = e" 2 «. 

In the remainder of this paper we will assume that the GB1 equation of state is valid throughout the 
star. Since Z is bounded it then follows from (^) that T and X both vanish at the center of the star. Using 
the translational freedom in R we may set R c — 0. The condition T(0) = then implies R_ = for £ = 
and that is a multiple of 7r for £ > 0. Moreover, to ensure that T is positive as R increases away from 

the center c<j_i?_ must actually be an even multiple of ir. Similarly, the condition X(0) — implies that 
oj+R+ must also be an even multiple of n. As already stated in M we may then set R± = without loss of 
generality. 

Let us pause for a moment to digest what is actually implied here. We have the curious situation that 
the mere assumption that the model can be continued to the center of the star implies the vanishing of 
one of the essential integration constants (D = R + — i?_ = 0). Technically this has happened because 
the function W(R) which gives Schwarzschild's radial variable is built up by two non-negative terms, r = 
W(R) — T(R) + X(R), which must both vanish separately at R = 0. Normally one would expect to use that 
integration constant to ensure the elementary flatness of the model at the center. So what about elementary 
flatness? 

It follows that the model in fact satisfies the elementary flatness condition |J. It seems that the GB1 
equation of state cannot support a conical singularity at the center. This result sparks off a number of 
questions. For example, what equations of state have this property of not allowing conical singularities? 
Also, since a conical singularity is not allowed, in what way does the model break down if we let say R + 
be nonzero? The answer to the latter question is that the equation of state itself becomes unphysical in 
the sense that u — > oo at some value of the Schwarzschild radial variable, r = > 0. At still smaller 
radii, r < r^, the matter becomes mechanically unstable, dp/ dp < 0. Thus by allowing a nonzero R + what 
happens is that the matter becomes unphysical before any breakdown in the geometry as we move towards 
lower r- values. 
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The radius of the star is defined as r s = W(R S )- The star mass can be computed by the formula jjj 

M= 1 1 [W(l~N-^ h J] s . (32) 

This leads to the following expressions for the physical parameters M, (3 and z s 

M = \r s [1 - Z(W, R ) 2 ] s , (3 = \ [1 - Z(W, R ) 2 ] s , 1 + z s = Z^^r)- 1 . (33) 

To evaluate these expressions at the surface we need to know R s . Using (|J) and (||) we find that R s = q/u) + 
where q is defined to be the smallest positive root of the equation sin(^x) = 77 sin a: where we have defined 
i] := tanh£. Thus q satisfies 

sin(£q) = 77 sing . (34) 

For Buchdahl's original model (£ = 0) we have q = ir and R s — 7r/(2 v / Ka). For £ 7^ equation ([mJ) cannot 
be solved explicitly. Now what about the size of 77 compared to £? In fact e -2 ^ = Z c is constrained by the 
requirement Z M < Z c < Z s . This leads to the inequality ( cf. ) 



Expressing this condition in terms of 77 leads to 



3/5 < = 35/5 < e- 2 < < 5 ^1 = ^ < 1 . (35) 



4^2 

T 2 



0<t<v<Vco:= -r—^- < 1/4 . (36) 



The allowed values of £ and 77 are shown graphically in figure |j. 

To determine the surface value of the radial coordinate we must solve equation (|34|) numerically for 
each set of values (£,77) in the parameter space defined by (13(f). However, this can be avoided by observing 
that 77 is an increasing function of q for all £. Conversely, q is an increasing function of 77. Therefore, 
we can use q as a parameter in place of 77 to represent the sequence of stellar models. It follows that q 
ranges from the minimal value q = corresponding to 77 = £ (see (|36|)) to a certain maximal value q = q^ 
corresponding to 77 — 77^, = (1— 4£ 2 )/(4— £ 2 ). The value q^ can only be determined implicitly by the equation 
sin^goo) = 77^ sinqoo. In the limit of small £ we find that q^ will be close to ir. Defining A := ir — q we 
have the relation 

(l-4e 2 )sinA oc = (4-£ 2 )sin[^- , (37) 

for Aoo. This equation can be used to solve for A^ as a power series in £ with the result 

Aoo = 4tt£ - 16< 2 + tt(79 + 10vr 2 )£ 3 + C(£ 4 ) . (38) 

An issue here is the radius of convergence of this series. At this point we will content ourselves with the 
observation that for small enough £ the above series solution ( |38| ) gives a good approximate solution to ( |37j ) . 
This is an indication that the radius of convergence is nonzero. For general £ it turns out to be more effective 
to solve ( |3"7j ) directly using a root finder routine. To obtain the value of A at the causal limit, A»(£), we 
proceed in an analogous way observing first that 77* = (1 — 5£ 2 )/(5 — £ 2 ), where 77* is defined by letting the 
sound speed at the center be equal to the speed of light, u c = 1. This leads to the equation 

(1 - 5£ 2 ) sin A* = (5 - £ 2 ) sin[£(7r - A,)] . (39) 

Again we solve this equation numerically using a root finding routine. The functions A oc (£) and A*(£) are 
plotted in figure |[ 

Although the parameter q is a good variable to describe sequences of models with £ 7^ it is not a good 
variable in the limit £ — * since q approaches the constant value n in that limit. We wish to obtain a variable 
which works well also in the £ — > limit. Since 77 is a good variable for the £ = models we look for a new 
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0.00 0.05 0.10 0.15 0.20 0.25 

Figure 4: The allowed £ — 77 region. The vertical dashed line represents the original Buchdahl model (£ = 0). 
The 77 = £ line corresponds to the low pressure Newtonian limit while the solid curve represents the most extreme 
relativistic models in the limit u c — > oo. From (|36|), the solid curve is given by the equation r\ = 7700 = (1— 4£ 2 )/(4— £ 2 ). 
The dotted line is the causal limit corresponding to models for which the speed of sound at the center equals the speed 
of light, u c = 1- Its equation is 77 = 77, := (1 — 5£ 2 )/(5 — £ 2 ). Models lying below the causal limit are entirely causal 
while models above this limit have an acausal region in the stellar interior. Models with £ > = 3 — 2v2 ~ 0.1716 
(the region to the right of the vertical solid line) are acausal throughout the stellar interior. For each given value of 
£ in the allowed range < f; < £co ~ 0.2087 there is a 1-parameter sequence of models represented by the vertical 
straight line which starts at the point (£,£) and ends on the solid curve. 




variable a which has the property a — > 77 as £ — > 0. To that end we observe from (|34|) that 77 = C'?/( 7r — l) 
when ( « 1. This motivates the definition a := £,q/(Tr — q). Then a — > 7/ in the limit £ — > as required. 
Also, (j is a nonnegative increasing function of q and hence is also an increasing function of 77. Finally we 
normalize a by denning a :— a/(Aa QO ) where (Too = £<7oo/ (it — <Zoo) = £(tt — A 00 )/A 00 . Then a — > 77 as £ — > 
and (Too = 1/4 for all £. We have 



* A °° (40) 



4( 7 r-g)(7r-A 00 ) ' 



and solving for g we obtain 



4 7 r( 7 r-A 00 )a 

A 00 +4(tt- A^ct v ' 



Expressing the causal limit in terms of the normalized parameter gives 

A 00 (7T-A >t ) 



4A*( 7 r-A 00 ) ' (42) 
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Figure 5: The function A 00 (^) is plotted in the interval < £ < £00 = (5 - v21)/2 « 0.2087. The function A*(f) 
is plotted in the interval 0<£<£,=3 — 2\/2 ~ 0.1716. The models in the range < £ < have a superluminal 
sound speed at the center for all values of rj. 



3.1 Expressions for radius, mass and surface redshift 

In the rest of the paper we set R± — 0. This implies that the model is regular throughout the interior of the 
star. To compute the radius we first use (ra) to write 



(£ v±~e 



W{R) 



2 VT^rf 



[f- 1 sin(w-R) + rj sm{bJ + R)] , for £ + 



(43) 



[u)+R + rj sm(w + i?)] 



for£ = . 



I 2 V^v 2 ' 

where £ := 1/^/ko, is a characteristic length scale. Evaluating this expression at the surface we find 



l(l + C)\A-£ 2 singsin(g g ) 
2£^sin 2 g -sin 2 (£< ? ) 
k£ 



for £ ^ 



for £ = 



(44) 



The radius is smallest in the limit 77 — * £ (or equivalently q — > if £ 7^ 0). It is evident from (|44|) that 
for BuchdahPs model the minimum radius is nonzero and is given by (r s ) m j n = (tt/2)£. The models with 
£ > on the other hand all have zero minimum radius as expected from their liquid-like equation of state. 
A contour plot of radii for the GB1 family is shown in figure || 
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Radius contours 
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5 

Figure 6: Contours of stellar radii are shown for the GB1 family. The radii are indicated in units of the characteristic 
length scale I. The solid contours are equidistant. The dotted contour has the value r s /£ = 7r/2 « 1.57. This is 
the minimum radius for Buchdahl's gaseous models. The models with £ > all have zero minimum radius reflecting 
their liquid-like equation of state. The causal limit as defined by ( [I^ ) is indicated by the grey line. Models below 
that line are everywhere causal, u < c. 



Surface redshift contours 
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Figure 7: Surface redshift contours are shown for all the GB1 models. The causal limit is indicated by the grey 
line. Models below that line are everywhere causal, u < c. 
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To compute the redshift we begin by writing down an expression for (3. For that we need 

I (cosh£) cos(uj-R) + (sinh£) cos(u;+i?) for £ + 1 , 



Wm = 



y cosh (+ (sinh£) cos(o; + i?) 



for£ = , 



(45) 



leading to 



(cosh () cos(£g) + (sinh £) cos q 



/ sin[(l + flg] 
sin[(l - £)g] 



for £ f 
for i = 



(46) 



This gives 




Qsin[(l + Qg] ] 
Osin[(l-0«]J 



The expression for the redshift becomes 



l + z s 



flsiii[(l-flg] 
Osin[(l + 0?] 




for £_ ^ 
for £ = 

for ^ , 
for f = . 



(47) 



(48) 



Comparison with the Doppler relation between redshift and velocity shows that for Buchdahl's model 77 
is exactly equal to the corresponding velocity. The maximum equivalent velocity for Buchdahl's model 
is therefore v = 77 = 1/4 corresponding to a maximum redshift of z s = \/5/3 — 1 ~ 0.291. Restriction 
to causal models (v — rj < 1/5) gives a maximum redshift z s = -\/3/2 — 1 w 0.225. A contour plot of 
surface redshifts for the GB1 family is shown in figure 0. The maximum value of the surface redshift can be 
computed numerically to be z ~ 0.31 which occurs for £ w 0.032. The corresponding compactness parameter 
is a — 1/ [3 ~ 4.8. Therefore there are no ultracompact models in the GB1 family. Turning now to the mass 
we find that it can be written as 



f ^(l + Ov/^F 



M 



4£y sin 2 q - sin 2 (£g) 



singsin(^g) f (1 - j) sin[(l + flgj \ , 



(49) 



[2(l + r ? ) N /T^ 



for £ = 



A contour plot of masses for the GB1 family is shown in figure 

The central values of the pressure and energy density are given by 

4a(?7 2 - i 2 ) 



Pc 



(i-e)(i+v) 2 ' pc (i-£ 2 )(i+7?) 2 

The central values of the adiabatic index and sound speed are 



4a[7 ? (2-3r ? )+£ 2 (3-2r ? )] 



7c 



2(1 -^(77+ a 



2\2 



77 +e 



W -C 2 )[l -4?7- (4-Tj)^ 2 ] 
The surface density is given by 



1 - 4t/ - (4 - 77)^ 



Ps = 



8a£ 



(50) 



(51) 



(52) 
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Mass contours 




0.00 0.05 0.10 0.15 0.20 

Figure 8: Mass contours are shown for all the GB1 models. The masses are given in units of the characteristic 
length, I. The causal limit is indicated by the grey line. Models below that line are everywhere causal, u < c. 



3.2 Pressure and density profiles 



The pressure and density functions are given in (| 
the normalized radial variable R := R/ R s we find 



as functions of Z. Expressing this function in terms of 



Z(R) = 



i - m 

1 + 4>(R) 



where 



<t>(R) := X/T = < 



£ sin(£g) sm(qR) 
smqsm(^qR) 

i] sin(irR) 
irR 



(53) 



fore^O , 
for f = . 



(54) 



These formulas together with (||) give p and p as functions of R. We wish to plot the pressure and density 
versus Schwarzschild's radial variable which was given in equation (^). The result for Buchdahl's model is 
shown in figure ||. A striking feature is the small variation in the radius. In going from a low mass Newtonian 
model with z = up to the most compact relativistic model with redshift z — 0.291 the variation in the 
radius is only about 3%. The profiles for a representative of the type II GB1 family are shown in figure Eol. 



4 Concluding remarks 

We have described the 2-parameter GB1 family of exact spherically symmetric relativistic stellar models 
with a physically reasonable equation of state. One of the parameters represents the scaling invariance 
while the second parameter can be interpreted as a measure of the stiffness. By restricting the range of the 
stiffness parameter, one obtains a subfamily for which the sound speed is subluminal throughout the stellar 
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r/£ rjt 

(a) Density profiles. (b) Pressure profiles. 

Figure 9: Density and pressure profiles for Buchdahl's model (£ = 0). The density and pressure are given in units 
of pa = K~^t~ 2 and po = k -1 ^ -2 . The profiles are given for four representative values of the surface redshift ranging 
from the mariginally relativistic, z = 0.01, to the most extreme relativistic CEIS6, Z — Zoo — s/5/3- 1 « 0.291. The 
radius takes values in the very narrow range, tt/2 < r/£ < 27r/\/l5, or approximately 1.571 < rjt < 1.622. 



interior. By varying the equation of state parameters it is possible to obtain a quite accurate fitting to the 
HW equation of state corresponding to the outer region of a neutron star. 

Being a realistic model of a neutron star near the surface the GB1 family is a good candidate for a static 
limit of a rotating neutron star model. The main difficulty in trying to spin up a static perfect fluid model is 
to preserve the perfect fluid property. Technically this means that the Ricci tensor should have the correct 



algebraic type (Segre characteristic Al[(lll),l)], see reference 22 ) and in addition contain only a single 



arbitrary function. A recent result |23f seems to indicate that the Newman-Janis algorithm would always 
destroy the perfect fluid property. This could mean that a rotating perfect fluid model would not have a 
Kerr exterior. 

The GB1 family is the only exact general physical SSS solution which is expressible in terms of elementary 
functions for all values of the integrations constants. For Schwarzschild's interior solution, for example, only 
the regular models can be described in terms of elementary functions while the singular cases involve elliptic 
functions. It should be noted that the singular models have considerable physical interest. First they can be 
used in the regions outside the core of the star. Also, they are needed for obtaining a better understanding 
of the space of solutions as well as of the possibilities for global behaviour. 
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(a) Density profiles. 

Figure 10: Profiles for the type II GB1 models (£ = 0. 
similar to those of the Buchdahl model except that the c 
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(b) Pressure profiles. 

0006). Axes and units are as in figure The profiles are 
[ensity is nonzero at the surface. 



Appendix 
Useful formulas 

First some definitions: 



5 - 1 



uj± := ^j2(8±l)Ka , i:=l/y/m, £ := uj- /uj+ = J - — - . (55) 

V o + 1 

We then have the following identities 

* = ^f, * + 1 = I^' 5 - l -T^e^ S-VP=T=^, (56) 

and 



= \ly/T=? , u-- 1 = l^V^ • (57) 

?/:=tanhC, e 2C = — ^ , coshC=^=L=, smb. C = — t^== • (58) 
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(a) Density profiles. 
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Figure 11: Profiles for type III GB1 models (£ = 0.06). Axes and units are as in figure 




(a) Density profiles. 



(b) Pressure profiles. 



Figure 12: Profiles for type IV GB1 models (£ = 0.2). Axes and units are as in figure ^. 
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